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We study the relaxation of a spin density injected into a two-dimensional electron system with 
generic spin-orbit interactions. Our model includes the Rashba as well as linear and cubic Dressel- 
haus terms. We explicitly derive a general spin-charge coupled diffusion equation. Spin diffusion 
is characterized by just two independent dimensionless parameters: 7^ and jd, which control the 
interplay between different spin-orbit couplings. The real-time representation of the diffuson ma- 
trix (Green's function of the diffusion equation) is evaluated analytically. The diffuson describes 
space-time dynamics of the injected spin distribution. We explicitly study two regimes: The first 
regime corresponds to negligible spin-charge coupling and is characterized by standard charge dif- 
fusion decoupled from the spin dynamics. It is shown that there exist several qualitatively different 
dynamic behaviors of the spin density, which correspond to various domains in the (7^,70) param- 
eter space. We discuss in detail a few interesting phenomena such as an enhancement of the spin 
relaxation times, real space oscillatory dynamics, and anisotropic transport. In the second regime, 
we include the effects of spin-charge coupling. It is shown that the spin-charge coupling leads to 
an enhancement of the effective charge diffusion coefficient. We also find that in the case of strong 
spin-charge coupling, the relaxation rates formally become complex and the spin/charge dynamics 
is characterized by real time oscillations. These effects are qualitatively similar to those observed in 
spin-grating experiments [Weber et al., Nature 437, 1330 (2005)]. 

PACS numbers: 

I. INTRODUCTION 

Spin-orbit coupled systems host an amazing variety of interesting effectSfi which are currently the subject of intense 
experimental and theoretical investigation. The interest in this field stems from the fact that spin-orbit coupling 
opens the possibility of controlling the electron spin^ using purely electrical means^i^i^i^ii with potential applications 
in spintronics^i^ and quantum computing^. The microscopic origin of the spin-orbit interaction in a low-dimensional 
system is the absence of inversion symmetry in the semiconductor crystal and confining potential. This leads to the 
appearance of the Rashba and Dresselhaus (linear and cubic) spin-orbit terms in the effective Hamiltonian.— li^ii^iii^ 
The spin-orbit interaction results in a variety of macroscopic effects such as antilocalization ; electric-field induced 
spin accumulation ; ^'''i spin-galvanic effects^^ etc. 

One particularly important aspect of spin transport in disordered systems is Dyakonov-Perel spin relaxationi^ 
The source of this effect is the precession of the electron spin around a momentum-dependent magnetic field. The 
electron momentum changes randomly due to electron scatterings off of impurities, which in turn randomizes spin 
precession. One should note here that a large spin relaxation rate is typically considered an undesirable effect from the 
experimental point of view. On the one hand, the Dyakonov-Perel spin relaxation rate is proportional to the strength 
of spin-orbit interactions. On the other hand, the existence of relatively strong spin orbit-interaction is needed for 
the appearance of non-trivial spin transport. It is therefore important to understand the interplay between different 
phenomena in spin transport and search for possibilities to reduce the negative effects of spin relaxation. 

We should mention here that the issue of spin relaxation has been extensively studied in experiments. In particular, 
Weber et al. have used the transient spin grating method to study spin dynamics in two-dimensional semicon- 
ductor systems i^^i^ii^. This experimental technique is based on optical injection of spin polarized electrons in a 
two-dimensional quantum well, which generates a spin density wave with a wave-vector in the plane of the two- 
dimensional system and the spin polarization perpendicular to the plane. Particularly interesting results revealed by 
these experiments^! were an enhancement of the charge diffusion coefficient and an unusual time-dependence of the 
spin grating amplitude, which exhibited oscillatory features in addition to the conventional exponential decay. 

Motivated by the aforementioned experimental results, we study theoretically spin transport and spin relaxation in 
a two-dimensional electron system with generic spin orbit interaction. We focus our attention on the interplay between 
different types of spin-orbit couplings and the effects of the spin-charge coupling on the spin relaxation dynamics. 
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We work within the framework of the diffusion approximation, which assumes that the spin relaxation length is large 
compared to the mean-free path. We find that depending on the relative values of the Rashba, linear and cubic 
Dresselhaus couplings one can have different dynamic spin-relaxation regimes and phenomena such as oscillations of 
the spin density in real space, anisotropic spin transport, and enhancement of spin relaxation times. We also study in 
detail the effect of spin-charge coupling on the diffusive dynamics. We find that it always leads to a renormalization of 
the effective charge diffusion coefficient, which gets enhanced compared to the "bare" diffusion constant. We also find 
that a strong spin-charge coupling may formally lead to a complex spin-relaxation rate and real-time oscillations in the 
spin and charge diffusion behavior, which arc qualitatively similar to the oscillatory behavior observed in experiment. 
We should note here that the real-time oscillations appear only beyond the formal regime of validity of the diffusion 
approximation. However, they give a tentative indication that the experimentally observed oscillations may originate 
from the spin-charge coupling. 

Our paper is structured as follows: 

In Sec. ini we introduce our model Hamiltonian for the non-interacting disordered two-dimensional electron system 
with spin-orbit interaction and derive a general spin-charge coupled diffusion equation. We introduce natural dimen- 
sionless units of length and time and show that in these units, the spin part of the diffusion equation is characterized 
by two dimensionless parameters, and 7^1. These parameters are constrained to lie within a circle of radius two 
in the (7^,7^) parameter space. We also introduce the matrix Green's function of the diffusion equation (diffuson), 
whose real time representation describes the spin relaxation dynamics of interest. 

In Sec. IIII[ we consider the spin sector of the diffusion equation (which is justified if the spin-charge coupling is 
relatively small). We focus our attention on the diagonal z-component of the diffuson, which describes the dynamics 
of injected out-of-plane spin density. We show that there exist different domains in the (7^^, 7£))-parameter space, 
which correspond to qualitatively different dynamic regimes. We explicitly derive the real space-time behavior of the 
diagonal component of the diffuson matrix in various regimes. We also discuss the spin dynamics as a function of 
momentum and derive the spin relaxation spectrum as a function of spin-orbit coupling parameters. In particular, we 
study in detail the "symmetry point" in the parameter space (7^ = 7-d), when the relaxation time becomes infinite 
in certain directions. We study the behavior of the spin relaxation spectrum in the vicinity of the symmetry point. 

In Sec. IIVI we discuss the effects of spin charge coupling on the spin relaxation spectrum. We show that the spin 
charge coupling may formally lead to a complex spin relaxation rate, which physically translates into an oscillatory 
behavior of the spin density. We discuss the validity of the diffusion equation approximation and the implication of 
our results for understanding spin transport experiments. Our conclusions and a discussion of several open questions 
are presented in the last Section [V] 



II. THE DIFFUSION EQUATION 



In this section we introduce the model and derive a general form of the spin-charge coupled diffusion equation. The 
main results of this section are equations pip and (|f 3p representing the most general form of the diffusion equation in 
an electron system with arbitrary spin-orbit interaction. In this and subsequent Sections, we use the units h ~ = \. 

We start with the following Hamiltonian, which describes the conduction band electrons in a III-V type semicon- 
ductor quantum well grown in the [001] direction (set as the z axis), 

i?=|^+h(p).a, (1) 
2m 

where m is the effective electron mass, = (ctx, o'y, o'z) is the Pauli matrix, and h(p) = (h^, hy, hz) are functions of 
the two-dimensional momentum p describing the spin-orbit interaction. We assume that the nonintcracting electrons 
are in the presence of a short-range impurity potential and we investigate the role of the spin-orbit interaction in the 
coupled spin and charge transport using the diffusion approximation. In general, the second term of the Hamiltonian 
([I]), which describes the spin-orbit interaction, contains both Rashba and Dresselhaus contributions. The Rashba 
interactioniiii^, which is due to the inversion asymmetry of the quantum well confining potential, has the form 

h'^ip) ^avF{Py,-Px), (2) 

where vp is the Fermi velocity, and a is a dimensionless coupling constant. Throughout this work we will assume 
that the spin-orbit coupling energy is much smaller than the Fermi energy, i.e. |a| ^ I. Consequently, we will neglect 
the change of the Fermi surface due to the presence of spin-orbit interaction. The Dresselhaus spin-orbit couplingfi^ 
arising from the lack of inversion symmetry of the semiconductor crystal, contains terms both linear and cubic in p, 

h^HP) = PlVF{p.,-Py), (3) 

h^^p) = ^AP,^{p^pI-p^pI), (4) 
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where /3i and P3 are dimensionless coupling constants with \(3j\ 1, and pf = mvp is the Fermi momentum. 

In the presence of disorder, the complete description of the coupled spin and charge diffusive transport at long 
lengthscales compared to the mean- free path is given by a set of linear differential equations that, in its most general 
form, can be written as 

(^^ - vv^^ = (-r^' + p^^"v™ + c^^ V) p, (5) 

where d/dt represents the time derivative, po is the charge density, while pi — px, P2 — Py, and pa ~ pz are the 
densities of the corresponding spin components. In equation ([5|) V represents the diffusion constant and can be 
expressed in terms of the mean scattering time r as I? = TVp/2, T^^ = {1/Ts)ij describes the Dyakonov-Perel spin 
relaxation's, py™ jg responsible for the precession of the inhomogeneous spin polarization, and the last term, having 
non-zero elements = C°', describes the spin-orbit mixing of spin and charge degrees of freedom. For convenience, 
throughout this article we will express times in units of spin relaxation time, Ts, and lengths in units spin relaxation 
length, Ls, unless otherwise stated. The quantities Tg and Lg depend explicitly on the parameters of the 2D electron 
system, 

^ 2fc^.[a2 + (/3i_/33)2+/j2]i/2- (6) 



2t 



(7) 



where g = 2vFkFT is a coefficient proportional to the dimensionless conductance. Using the units (j6|7p , we can write 
the coefficients from Eq. ([5]) as dimensionless coupling constants. For example, the diagonal contributions to the spin 
relaxation terms are T^x ~ ^yy = 1 and T^z ~ 2. The final matter of convenience concerns the paramctrization of 
the spin-orbit interaction. This interaction can be described in terms of the original coupling constants (a, /3i, P^) or, 
alternatively, in terms of a new set of parameters, (r,7/j,7£)), defined by 



a'+l3',+ Pi (8) 



of /33. We derive the diffusion equations using the standard density matrix formalism's, used previously in Refs. 



2a 2(/3i - /?3) 

The parameter T represents a measure of the overall strength of the spin-orbit interaction, while and characterize 
the interplay between the Rashba and Dresselhaus contributions. By convention the sign of F is the same as the sign 

daiillillOl. Eq. © can by re- written in the form 

- %] p, = 0, (9) 

where, for a homogeneous system, Sij — 11^ (0;, k) = {—iuj + k^)Sij + F*-' -I- P^^"^km + C^^k and summations over 
repeated indices are assumed. The matrix II represents the kernel of the diffusion equation and can be expressed in 
terms of retarded and advanced Green functions as 

%(^,k) = / T^Trfa, G«(u;/2, k + q/2) a, G^(c./2, k - q/2)| , (10) 

where G^^'^^ is the 2x2 retarded (advanced) Green function matrix averaged over disorder, ct^ are the Pauli matrices, 
vp is the density of states at the Fermi energy and Tr{. . .} involves the trace over spin indices. The diffusion 
approximation represents the low frequency, long wavelength limit of Eq. (jlOp . This approximation is valid as long 
as the scattering time is much shorter than the spin precession time, i.e. in the weak spin-orbit coupling limit 
TvphpT <C 1. Taking the appropriate small frequency and small momentum limit, we obtain 



i - n(tJ,k) 



f s-1 - ig{pokx + fJ-iky) ig{pikx + poky) \ 

-ig{pokx + piky) s \ir1d - i[lRkx + loky) 

igipikx + Poky) \lRlD s - ii^okx + jRky) 

\ iilRkx + Joky) iilnkx +lRky) s + 1 ) 



(11) 



where 1 is the 4x4 unit matrix, s = —iuo + fc + 1 and the spin-charge coupling is characterized by the parameters 

{?>P^^P,)(a^^Pl+Pl)-P,Pl 

1^0 = ■ 



Va2 + (/3i-/?3)2 + /3|- 



^ (12) 
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In addition, it will be convenient to consider the diffusion problem within a coordinate system that is rotated in 
the x-y plane counterclockwise with 7r/2 relative to the initial coordinate system^. The rotated coordinates are 
r_|_ = (y + x)l\f2 and r_ = — a;)/V2, while the new components of the spin density will be p+ = (pj, + px)l\f^ 
and p- = {py — px)/V^, in addition to pz and the charge density po that remain unchanged. Relative to the rotated 
coordinate system, diffusion is described by the operator 



i-n('')(tj,k) 



i^{pQ+Pi)k+ s-^jR-jD 



\ 

-«(7-R + 7D)fc+ 

-«(7-R^ - 7D)fc- 
s + 1 j 



(13) 



where the four columns correspond to po, /0+, p~^ and p,, respectively. Equations (|lip and p3p are the main results 
of this section. 

The physical problem that we will be solving on the basis of Eqs. (fTTl[T3l) is to describe the transport of an injected 
spin density in the presence of a general type of spin-orbit coupling. Assuming that a spin (or charge) density with 
components pi(0, r) was injected into the system at i = 0, the density profile at times t > will be 



P»(<,r) 



y" dv' Aj (<,r,r')pj(0,r'), 



(14) 



where D is the diffuson, i.e., the Green's function of the diffusion equation. The physical meaning of the diffuson 
can be understood if we assume that initially we inject a i5-like, density profile pj(0,r) = <5(r). The spin and charge 
densities at any later time are given by the matrix elements of the diffuson, pi{t,r) = D,y (t, r). For a homogeneous 
system the diffuson is the inverse of the diffusion kernel. 



D 



1 - n 



(15) 



Consequently, solving the transport problem for an injected spin density implies inverting the matrix in Eq. (jlip or 
(|13p and performing the appropriate Fourier transforms. We notice that the spin and charge dynamics is controlled by 
the poles of D{lo, k). These poles determine four relaxation modes iLUa(\i) = l/TQ(k), which follow from the equation: 



A(w,k) = det[i - fl(s,k)] = 
A time dependent matrix element of the diffuson can be written in general as 
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D,y(t,k)=^A,(k) 



(16) 



(17) 



1=0 



where the amplitudes Ai{\i) are functions of momentum. We should emphasize here that the long time relaxation of 
the injected spin density is uniquely determined by the lowest minimum of the relaxation rate modes. In a standard 
diffusion problem, the charge relaxation rate is l/ro(k) = Vk^ (where T) is the diffusion constant, which is equal to 
one in our special units) and the minimum obviously corresponds to k = 0. By contrast, in the presence of spin-orbit 
interaction, the k-dependence of the relaxation rate is more complicated, with possible minima at finite momenta. 
The detailed dependence of the relaxation rate on the spin-orbit interaction is analyzed in the following sections. 



III. THE SOLUTION OF THE TIME-DEPENDENT DIFFUSION EQUATION: THE SPIN SECTOR 

In this section we discuss the spin relaxation spectrum and derive the diffuson in momentum space IIII Al and real 
space HITB] The main results of section UlI Al arc Eqs. (PTI describing the spin relaxation spectrum and Fig. [U 
which summarizes the properties of the spectrum. The behavior in the vicinity of the symmetry points characterized 
by an infinite spin relaxation time is illustrated in Fig. [2l Eqs. (|30p and (|33ll34p describe spin diffusion in real space. 



A. Momentum space picture 



The diffusion approximation leading to equation (jlip is rigourously valid only in the weak spin-orbit coupling regime 
characterized by TvpkpT <^ 1. In this limit, the spin-charge coupling is small, 5|/^i.2| ^ 1, and can be neglected 
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in the leading approximation. The goal of this section is to determine the role played by the interference between 
different spin-orbit interaction terms in spin diffusion. In particular, motivated by experiments that probe spin 
dynamics optically, we concentrate on the properties of the out-of-plane component of the spin density. To address 
this problem, we have to determine the element D^z = [1 — n~^]zz of the Green's fimction for the diffusion equation 
(the diffuson). The general expression for the diagonal z-component of the diffuson in the absence of spin-charge 
coupling is 

i?.,(s,k) = [s2_(^^^^/2)2]/A(5,k), 
where A(s, k) is a third order polynomial in s. 



A(s,k) = s^{s + l)-s 



222 / 2 2 \ ( '~iRlD \ ■' 

+ IrId^ + 1R.1D [1r + Id) kxky - ( 1 



(18) 



The spin dynamics is determined by the three spin relaxation rate modes iujji^L) ~ — (k) obtained by solving 

the equation A(s,k) = 0. For later reference, let us identify the modes that control the spin dynamics, together with 
the charge mode, according to their zero momentum values. 



zcjo(O) 
iwi(O) 

IW2(0) 



= 0, 

= 1 

= 1 

- 2, 



IRID 

2 ' 

IRID 



(19) 



where the mode iojq is responsible for charge diffusion and has the standard form ia-'o(k) = fc^ if we neglect spin-charge 
coupling. Let us first focus on the experimentally relevant special case of momentum parallel to the [110] or [110] 
directions corresponding to k±. Inverting the matrix from Eq. (|13[) we obtain. 



Dzz{s,k±) 



s± 



Trip 
2 



S 1± 



2 , 



4 



ill 



±7d)2±^ 



(20) 



Notice that, for k-vectors oriented along these special directions, one of the modes icui or iuj2 does not contribute to 
the dynamics and, consequently, the problem simplifies significantly. The time dependence of Dzz can be obtained by 
Fourier transforming (|20p with respect to the frequency. This time dependence is uniquely determined by the roots 
of the secular equation A(s, k±) ~ 0, i.e. by 



iwi(fc_) = 2+^- 1 [ ^)~) +4fc-(7i? -7I))^ 



-kl 



4 2 

iRlD 1 



iuj2{k+) = - 

iuJ3{k±) = ^ + k±±^— + -]^[lT^^) +4fc4(7ii±7D) ■ 



(21) 



As mentioned above, the modes ia;i(fc+) = 1 + k'^ — Jr.Jd/'^ and iuj2{k-) = 1 + k'^ + 7/j7, 
the dynamics of the out-of-plane component of the spin. The explicit time dependence of D 
direction is 



Dzz{t, k+) = - [e-*"^('=+)* + e-'''^^('^+)* 



Ir.Id/'^ do not contribute to 
for k-vectors along the 

(22) 



A similar expression can be obtained for Dzzit,k-) by replacing fc+ with fc_, 1^2 with iuji and changing the sign of 
7£). The spin dynamics at large times, t 3> t^, is controlled by the minima of the relaxation rate modes iujj{k±). 
The mode itu^ has always the minimum at fc = 0, 1^2(0) ~ 2. On the other hand, the position and the values of 
the minima for the other two modes depend on the spin-orbit coupling parameters. The mode iwi has a minimum 
iwi(O) = 1 — 7_r7d/2 at zero wave- vector if 7^ + Jd ~ ^IrId/^ < 1. Otherwise, the dispersion curve has a local 
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FIG. 1: Diagram summarizing the the types of minima characterizing the lowest energy relaxation rate modes iuji and iix)2 
(see Eq. ((21}). The spin sector is uniquely described by the dimensionless spin-orbit coupling constants 7_r and 70. These 
parameters are proportional to the Rashba interaction and to the difference between the linear and cubic Dresselhaus couplings, 
respectively, and they satisfy the physical constraint 7^ + 7|j < 4 (see Eq. (|8])). In the red zone, both modes icoi and iLJ2 
have the minimum at k = 0. In the blue (yellow) regions, the minimum of the iwi (1102) mode is at zero momentum, while 
the minimum of the ii02 (ii^i) mode is at a finite wave- vector (k'L) given by Eq. (|23|l . In the green zones, both modes have 
a finite k- vector minimum at (for iuji) and (for 10^2). The points R and R' correspond to the the pure Rashba case 
(7«, 7d) = (±2, 0), while the segment DD' corresponds to the pure Dresselhaus spin-orbit interaction, jr = 0, 70 £ [—2,2]. 
At the end points 70 = ±2 the cubic Dresselhaus term vanishes. The spin dynamics is isotropic for 7^70 = 0, i.e. along the 
RR' and D'D segments. Sj and S'j are special symmetry points characterized by a zero minimum relaxation rate iujj{k'^), i.e. 
an infinite spin relaxation time. 



maximum at A: = and the minimum occurs at a finite wave- vector . Similarly, the mode iuj2 has cither a zero wave- 
vector minimum, if 7^ -|- 7!, -|- 57fl7D/2 < 1, or otherwise a finite k- vector minimum at fc^. The finite wave- vectors 
that describe the positions of the minima are 



i4i^ v^'""!:::"r^'' . (23) 



T in. I u ^ ^ 
2 

2|7fl ± 7£) 

while the frequencies corresponding to the minima in the dispersion curves are 

3 IRID (1 + ^)' ilR^lD? 



2 4 4(7^-715)2 4 



. ^ + ™ - - . (24) 

^ +^ 2 4 ^{^R + iD? 4 ^ ' 

A summary of the conditions necessary for the existence of these minima in terms of the spin-orbit coupling parameters 
7_R and 7d is presented in Fig. [1] In the red domain the condition 7^ -|- 7|j ± 57^72:1/2 < 1 is satisfied and all the 
dispersion curves have minima at zero momentum. The blue regions are characterized by a finite wave- vector minimum 
of the iL02 mode along the kj^. direction and a zero momentum minimum of the iuji mode, while the yellow sectors 
correspond to the [iuii{k^_), *wi(0)] pair of minima. Finally, in the green regions both modes, iuji and iuj2 have finite 
k- vector minima corresponding to k'^_ and k'i^, respectively. Notice that spin diffusion is isotropic only if 7j?7d = 0, 
i.e. for pure Rashba spin-orbit interaction (/3i = /^a = or (7_r,7_d) = (±2,0), the points R and R' in Fig. [T]), 
pure Dresselhaus spin-orbit coupling (a = or 7^^ = 0, the segment D'D in Fig. [T]), or a combination of Rasba and 
Dresselhaus contributions with the linear and cubic Dresselhaus terms being equal {(3i = (3^ or 7/3 =0 and a 7^ 0, 
the segment RR' without the end points). In all the other cases the spin diffuses anisotropically as a result of the 
interference between the spin-orbit coupling terms. 

To be more specific, let us consider now the case > 0, > 0, i.e. the quarter defined by the points R, O, 
and D in Fig. [T] Wc assume that a uniform out-of-plane spin density was injected into the system and ask about its 
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time evolution. The answer can be obtained using Eq. ((22)) for fc+ = 0. For a uniform system, the out-of-plane spin 
diffusion equation that takes the trivial form ^^^(ijO) = exp(— 2t), which corresponds to ia;3(0) = 2 representing the 
uniform spin relaxation rate of the S'^-component. If we consider now the limit of small non-zero momenta an expand 
all the quantities up to second order terms in k, we obtain 

D,,{t,k) « e ^C-)* + + fc^e— ^(•^)* + ^C-)*, (25) 

(1 + ^) 

where the dispersion relations at small k-vcctors take the form 

zc.2(k) « + - ki, (26) 



1 + 


7H7D 

2 


ilR 


+ 1d)^ 


1 - 


IRJD 


{jR 




1 + 


IRJD 

2 



lW3(k) W «W3(0) + fc + ^ , fc_ + —— k 



1 - 



2 



where /c^ = fc^ + fc^ = fc5]_ + fc^ . The initial decay of the spin density is dominated by the third term in Eq. (|25p that 
has the form exp[— (2 + I?+fc^ + 'D^K^)t]. Here, 'D± are effective diffusion constants along the k± directions. This 
result shows explicitly that the diffusion of the out-of-plane component of the spin is, in general, anisotropic, and that 
the anisotropy increases rapidly as we approach the symmetry point Si . Another consequence of these relations is 
that the asymptotic behavior at t 3> is controlled by the iuji mode, unlike the uniform case described solely by iuj^. 
In the isotropic case, the relaxation times corresponding to these modes differ by a factor of two, Ti/t^ = iuj^/iuji = 2. 
As we approach the symmetry point, this ratio diverges, ti/ts = 4/(2 — 7_r7d)- At the same time, the pre-factor 
of the first term in Eq. (|25p vanishes as we approach the OSi line (sec Fig. [T|) characterized by equal Rashba and 
Dresselhaus couplings, or if the wave- vector is oriented along the fc_|_ direction. Wc conclude that in general the 
asymptotic spin relaxation at small (but not- vanishing) k- vectors differs strongly from the spin relaxation at A: = 0. 
In addition, the small wave- vector behavior is, in general, anisotropic as a result of the interference between different 
spin-orbit coupling terms. 

Next, we focus on the finite momentum minima of the dispersion relations ([23]) and ((24)) . The existence of these 
minima for the itoi and iuj2 modes implies an increase in the spin lifetime relative to fc = for the corresponding 
mode by a factor iuji(^2){^) /^'^i{2){k^) ■ For example, in the pure Rashba case wc obtain k'^ = \/T5/4 in units of l/Lg, 
or kj. ~ akpVi^/'^ in standard units, and an increase in the lifetime by a factor of 16/7, in agreement with previous 
calculation^Zi^ based on the Rashba model. 

Another important special case is represented by the symmetry points characterized by equal Rashba and Dres- 
selhaus couplings, more precisely defined by = |7£i| ~ (the Sj points in Fig. [T]). At a symmetry point 
the minimum in the dispersion curve vanishes, iuj2{i){V2) ~ 0, reflecting the absence of spin relaxation^^i^ at the 
corresponding wave-vector, fc*J_ ~ \/2 (or /c° = -\/2 if 7fl7D < 0). The existence of a finite k-vector minimum in the dis- 
persion curves is illustrated in Fig. [2| The black curve with circles corresponds to the symmetry point 7^^ — ~ 
and has a gaplcss minimum at k^ = ^/2/ Lg- As we move away from the symmetry point, a gap develops (red and 
green lines). How is the gap behaving as we vary the coupling parameters in the vicinity of the symmetry point? To 
answer this question, we define d-^ as the distance in the parameter space {'^b.tId) from the symmetry point 

ci^ = [(7fl- V2)2 + (7^-V2)2]i/2. (27) 

If we move away from the symmetry point by changing the ratio between the Rashba and linear Dresselhaus interac- 
tions in the absence of a cubic Dressselhaus term, i.e. by varying the ratio a/ (3i and keeping (5^ = (see the inset of 
Fig. [21 red curve), the gap remains very small in comparison with the zero momentum frequency iL02{Q) « 2/r5. In 
contrast, the gap develops rapidly if we move along a direction in parameter space corresponding to = Id (green 
line). This observation suggests that, in order to observe experimentally the "absence" of spin relaxation at a certain 
finite momentum, the minimization of the cubic Dresselhaus interaction should be a major concern. In addition, we 
have to remember about the small momenta relaxation times given by equations ()25p and ()26p , in particular the one 
associated with the iuji mode. To summarize the results concerning relaxation times, we refer to the diagram in Fig. 
[TJ For a uniform system (k=0), the only mode that contributes to the dynamics is zws, and the corresponding spin 
relaxation time of the out-of-plane component Sz is Ts/2, regardless of the spin-orbit coupling parameters. The system 
can have a slower relaxation at non-zero momenta that can be either arbitrarily small or of order l/Lg, depending 
on the values of the coupling parameters 7/j and ^jj. In the red region, all three modes iwj(k) have the minimum at 



8 




FIG. 2: Momentum dependence of the relaxation rate mode ii02{k+) for spin-orbit coupling parameters in the vicinity of the 
symmetry point '^b. = 'yo = \/2 (point S\ in Fig. [T]). In this regime, the iuj2 mode has a finite wave- vector minimum along 
the fc+ direction. At the symmetry point, the dispersion curve becomes gapless (black curve and circles), meaning that the 
spin relaxation time at k% — \/2/Ls is infinite. As we depart from the symmetry point, a finite gap opens. The spin-orbit 
coupling parameters for the red curve correspond to a point on the RSiD arc in Fig. [1] while the green curve is for a point 
on the SiO segment, both points being in the vicinity of Si. The gap size as a function of the distance from the symmetry 
point [see Eq. (|27|l for the definition of dry] is shown in the inset. The red curve correspond to a spin orbit interaction with 
vanishing cubic Dresselhaus contribution /Ss = (the RSiD arc), while the green curve was obtained for 7h = 70 < a/2 (the 
SiO segment). To observe a substantial enhancement of the spin relaxation time at finite momentum, i.e. to have a small gap 
in the relaxation rate curve, the cubic Dresselhaus spin-orbit interaction /Js should be minimal. 



k = and the largest spin relaxation time, Tmax{0) = Ts/{1 — |7ij7£)|/2), is obtained at small k-vectors. In the blue 
(yellow) region, in addition to this small momentum minimum, we have the minima at fc*^ (k^) given by equations (|23p 
and (j24[) and to determine the slowest spin relaxation we have to compare T„iax (0) with T„iax 

Finally, in the green regions, the largest spin relaxation time is associated with one of the minima at k^.. The diffusion 
is isotropic if 7_r7d = and anisotropic otherwise, especially in the vicinity of the symmetry points. 



B. Real space picture 



We turn now our attention to the real space behavior of D^z, corresponding to the diffusion of a delta function- like 
spin density injected in the origin at < = 0. While the analytical treatment of the general case its rather complicated, 
we can easily solve special cases, like the isotropic case {'^r^d = 0), or the symmetry point 77^ = ±713 = \/2. 
In particular, for the isotropic case we obtain from Eq. ()lip . after inverting the matrix and performing a Fourier 
transform with respect to the frequency. 



^..(i,k) 



^1 + 472^2 



^1 + 472fc2 



-(f+fc=-i^l+472fc2^( 



(28) 



where 7 = 7/j if — and 7 = 7^) if we have a pure Dresselhaus spin-orbit coupling. The real space dependence is 
obtained after performing another Fourier transform, with respect to momentum, which reduces to the one dimensional 
integral 



1 

47r 



dk k Jo(fcr)D,,(t,k), 



(29) 



where Jo(fcr) represents a Bessel function of the first kind. The time dependence of Dzz{t,v) at distance r form 
the origin is shown in Fig. [3] for the particular case of pure Rashba coupling 7 = 2. The curves represent the time 
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t 

FIG. 3: Relaxation of tlie out-of-plane spin density in the pure Rashba case. The spin density p,, having an initial 5-like profile 
in real space, is injected in the origin at i = 0. The curves show the time evolution pz{t,r) — Dzz{t,r) at various distances r 
from the origin. The large time behavior is described by Eq. (|30[l for 7 = 2. Time is expressed in units of Ts and length in 
units of Ls. 



evolution of the out-of-plane spin density at different positions in the plane after spin was injected in the origin at 
t = 0. Notice the oscillations in the direction of the spin polarization (see Fig. 2]). The analytical expressions for these 
oscillations can be derived in the large time limit using a saddle point approximation for the integral over momenta 
[more precisely in the limit I/ts ^ [r/Ls)^]- For the isotropic case this asymptotic behavior has the form 



2 4 
- t 



where 7 = 7_r or 7 = 7^1, as appropriate. The asymptotic behavior is qualitatively different for |7ij| > 1 (or \'^d \ > 1), 
the green regions in Fig. [1] and |7fl| < 1 (or < 1), red zone. In the first case, Dzz{t,r) oscillates in space 
at large distances away from the origin, r/Ls ^ 1, with a period A = 2i:/{Lsk^) and an amplitude that decays 
as exp[— iix'(fc°)t]/Vrt, where A:" = fc'J^ = fc" and iuo{kF') ~ iuJi{k^) ~ iuJ2{k^) are the parameters for the finite 
momentum minima of the dispersion relations described by equations (|23p and (|24p . By contrast, in the red region, 
spin diffusion at large times, I/ts 3> (r/Lg)^, is independent of position. In addition, the spin relaxation time has 
the value r,, independent of the spin-orbit coupling, and the time dependence of the pre- factor is l/t"^, instead of the 
1/Vi dependence for the oscillating case. 

For "fR^D 7^ the dynamics of the Sz spin component is anisotropic. The simplest analytical treatment of such 
a case obtains at the symmetry points |7/j(d)| = V^- Starting with Eq. (|13p . we get, after inverting the matrix and 
performing a Fourier transform with respect to the frequency, 



D 



ZZ 



{t, k) = ie-[(v^+'^+)'+'^-]* + }_^-[{V2-k+ f+kl]t^ (^3-^^ 



2 2 



where we assumed that ^b.1d = 2, which corresponds to the symmetry points Si and S'l in Fig. [T] A similar equation 
can be obtained for ^r^d = ~2 by exchanging fc+ and fc_. This equation does not acquire any spin-charge coupling 
corrections, as the coupling vanishes at the symmetry points, /^o = Mi = 0, regardless of the strength of the spin-orbit 
interaction. Next, after performing the Fourier transform with respect to momentum, we obtain the exact real space 
dependence at the symmetry point, 

DUt, r) = ^ e-^ cos(V2 r+), (32) 
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FIG. 4: Real space dependence of the out-of-plane spin density in a system with isotropic spin-orbit interaction. The curves 
represent Dzz{t,r) muhiphed by the exponential factor exp(f/rm), where t = 36 (in units of Ts) and = (672 — 74— l)/(472) 
represents the maximum relaxation time corresponding to the finite k'^ minimum of the relaxation rate. At short distances, 
r <^ \/t — 6, the solutions are described by Eq. (|30|l . The period of the spatial oscillations increases as we decrease 7: 7 = 2 
(pure Rashbe or pure linear Dresselhaus interaction) - red line and circles, 7 = 1.2 - green line and triangles, 7 = 1- blue line 
and stars. Notice that, as we approach the non-oscillatory regime characterized by 7 < 1, the spin density becomes independent 
of distance in the asymptotic limit t ^ (see Eq. (|30|) ) and the 7=1 curve). The inset shows the spatial dependence of the 
spin density of a pure Rashba system (7 = 2) at different times. Again, we multiplied the results by the exponential factors 
exp{t/rm)- Notice that in the limit described by Eq. pop the curves are proportional, with proportionality factors 1/3 (black 
curve and diamonds, t — 9), 1/4 (cyan with squares, t = 16) and 1/6 (red line and circles, t = 36). 



where = + and we assumed that jr'Jd ~ 2. Because of the gapless minimum of the iLU2{k+) mode, the 
large time decay is proportional to 1/t, rather than exponential. In addition, the spin density oscillates along the re- 
direction, due to the finite momentum minimum, while it becomes independent of in the large time limit. For 
a general anisotropic case, the asymptotic behavior is determined by the minimum of the lowest relaxation rate. A 
characterization of the locations of the lowest minimum in the plane of the spin-orbit coupling parameters is provided 
in Fig. [SI We determine the asymptotic spin dynamics using a quadratic approximation for the dispersion curves in 
the vicinity of the minima. The approximate solutions of the secular equation A(s, k) for wave- vectors with arbitrary 
and small fc_ components are given in Appendix|3 Assuming that jajo > 0, the zero momentum minimum of the 
ioji mode represents the lowest energy contribution for a system with spin-orbit coupling parameters corresponding 
to the pink zone in Fig. \5\ The asymptotic S^-spin diffusion in the zone I of the parameter space (see Fig. [5]) is 
described by 



where the coefficient Cq^'' is given by Eq. (jA4p . Eq. ([55)1 is valid for times larger than the spin relaxation time, 
t ^ Ts, when ^ Jd- When the couplings are equal, the iwi mode whose fc = minimum generates the asymptotic 
behavior described by Eq. (|33[) does not contribute to the spin dynamics. Instead, for = corresponding to the 
OP segment in Fig. [5l the lowest minimum is the zero momentum minimum of the 10^2 mode and it will generate 
an expression similar to Eq. (|33p . but with r± — > and "fu ^ ^1d- Notice that in this regime spin diffusion is 
anisotropic only for times of order or smaller and becomes isotropic at longer times. Finally, for the anisotropic 
regime controlled by a finite momentum minimum in the dispersion relations, the blue and yellow regions and the 
segment PSi in Fig. [51 the asymptotic expression for Dzz can be obtained using the quadratic expansions given in 
Appendix [^ in particular equations (jA7[) and (|A8|) . We obtain the asymptotic expression for the diffuson in the 



11 




FIG. 5: Diagram summarizing the asymptotic relaxation of the out-of-plane spin density in real space for a quarter of the 
parameter space (7_r > 0, 7/3 > 0). In the pink zone (except the segment 05*1 corresponding to 7h = 70) the solution is 
anisotropic at short times and becomes isotropic at large times. This solution is controlled by the zero momentum minimum 
iuji{0) = 1 — 7h7d/2 and has the asymptotic expression given by Eq. (|33|) . For the segment OP, the solution is controlled by 
the zero momentum minimum iuj2{0) = 1 + 7^713/2 and has the same general characteristics as the neighboring pink zone. In 
the blue region the solution is anisotropic at all timescales and features real space oscillations along the direction. The long 
time dynamics is controlled by the finite momentum minimum iuj2{k%) and the asymptotic behavior is described by Eq. (|34|l . 
The same type of solution is valid for the PSi segment. Finally, in the yellow zone, the solution has real space oscillations 
along the '-' direction and an asymptotic behavior controlled by the finite momentum minimum icjilk'L). Notice that, as we 
approach the axes 7^ = and 7d = the spin relaxation, as described by Eq. ((28}, becomes isotropic. We stress that the 
regimes represented here are well defined in the long time limit. At some intermediate timescales, and especially near the 
zone boundaries, the solution may exhibit mixed characteristics when the minima of two distinct relaxation rate modes have 
comparable contributions. 



regime described by zone II of the parameter space (see Fig. [S]) 



Airt 



^(2) 



r(2) ^(2) 



: exp 




ii02{kl)t 



cos{k^r^), 



(34) 



where Z?^?,^ = 1 — (1 — ^r^d/'^)/ {ir + 7d)^, B^jm and C^^o are given in Appendix El and iu;2(fc^{.) represents the 

k^ minimum of the iuj2 mode given explicitly by equations (j23p and (j24p . Equation p4p describes the large time 
dynamics of the out-of-plane spin component when the ^0^2(^4.) represents lowest energy minimum, i.e. for the blue 
zone and the segment PSi in Fig. [S] A similar expression can be obtained when the dynamics is controlled by the 
iuji{k'^) minimum (yellow zone III) using the correspondence rules 1002 — > ioJi, r_ and ^ ~1d'- 



D^f(i,r) 



D 



(2) 



r,(2) ^(2) 



; exp 




ioJi{k°_)t 



cos(fc" r_). 



(35) 



where D; 



fo' - 1 - (1 + lRlD/^)/hR - 1d)^ and the coefficients b'^^J and cj^^ are given in Appendix El Notice 

that Eq. (|34|35p reduce to Eq. ([3T|l at the symmetry points jr ~ — ±^2 where it becomes exact. Within this 
regime, the solution of the diffusion equation is characterized by real space oscillations with a period A = 2'KLs/k^ 
determined by the position of the minimum in k-space while, at times much larger than r^/L^, the solution becomes 
independent of r_. The relaxation time. Ts/(iuj2{k^)), can be much larger than the relaxation time in the regime 
controlled by the zero momentum minimum, Tg/ {iuii{Q)), and diverges as we approach the symmetry points. On 
the other hand, when we approach the 77^ = or 7D = axes, the solution becomes isotropic and some of the 
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approximations used in deriving equations (|33p and (j34p are not longer valid. Instead, the asymptotic behavior in the 
isotropic regime is described by Eq. ([50]) . 

IV. THE ROLE OF THE SPIN-CHARGE COUPLING 

In this section we study the effects of spin-charge coupling on spin and charge diffusion. The main results of 
this section are equations ([57)) and (|42p describing the behavior of the relaxation rate modes in the presence of the 
spin-charge coupling. 

In the absence of spin-charge coupling, the charge diffuses according to the standard diffusion equation, while a 
spin density characterized by a given k-vector decays, in general, as a sum of three exponentials. Recent spin grating 
experiments^! have shown that, in certain conditions, the time decay of an optically injected spin density may include 
an oscillatory component. In this section, we explore the possibility that the appearance of these oscillations is due 
to spin-charge coupled dynamics. To this end, we re-derive the Green's function of the diffusion equation starting 
from the polarizability matrix Eq. (jlip . or from Eq. (|13|) . but without neglecting the spin-charge coupling terms 
proportional to fiQ and /xi. Inverting the polarizability matrix produces, in general, a complicated expression for 
the diffuson. However, there are several special cases in which significant simplifications occur. These cases involve, 
on the one hand, special values of the spin-orbit coupling parameters, namely the pure Rasliba interaction, 7/5 = 0, 
7i?, = ±2, the pure Dresselhaus coupling, 7^, = 0, and the symmetry points = ±>/2, 7_d ~ ±a/2. On the other hand, 
the inversion of n(ti;,k) can can be performed analytically for arbitrary values of the spin-orbit coupling parameters 
along the special directions in momentum space corresponding to and Our strategy is to obtain explicit 
analytical expressions for these special cases, and use them, in conjunction with numerical results for arbitrary sets 
of parameters, to extract and characterize the main features of the general solution. We start with the isotropic case 
corresponding to a pure Rashba coupling, or a pure Dresselhaus spin-orbit interaction. The modes that control the 
coupled spin-charge dynamics are the solutions of the secular equation 

A{s, k) = [s{s + 1) - j^k'] [s{s - 1) + g'fifk'] = 0, (36) 

where A(s, k) is the determinant of the II(s, k) matrix, s = —iui -I- 1 -I- fc^ and 7 = ±2 for the Rashba case, or 7 = 7D € 

[—2, 2] for a pure Dresselhaus spin-orbit interaction. The spin-charge coupling is measured by the parameter and 

1/2 

we have ^/ = /ii = sign(Q;)Q!^ in the Rashba case and fJ-i = fJ-o ~ [(3/33 — Mifil — (3f) — /3i/3|] / [(/3i — /^s)^ + /?|] 
for the Dresselhaus coupling. The full expression of the diffusion matrix is given in Appendix [B] The real time 
behavior of the diffusion propagator can be obtained after performing a Fourier transform of (jBl[) with respect to the 
frequency. This time dependence is determined by the roots sj = —iujj + 1 -f A:^ of equation (j36p . and has the general 
form (fTT)) . For the pure cases described by the equations ([M)) and (|Bip . the four modes that control the coupled 
spin-charge dynamics are 



' iwo(k) = 


1 
2 






- ig'^fifk'^ 


icji(k) = 


1 

2 


hP + 




- ig'^fi'fk'^ 


iw2(k) = 


3 
2 


hP- 




f A'j^k^ 


^ iw3(k) = 


3 
2 






f A'y^k'^ 



where g = 2vpkFT, fc^ = k^ + ky and the k-vector is expressed, as usual, in units of l/Lj. Equation ([57)) has several key 
features. First, we notice that the mode responsible for the charge dynamics in the absence of spin-charge coupling, 
namely iujQ, couples with only one spin mode, iwi, while the other spin modes, 11^2 and iuj^ remain unaffected by the 
spin-charge coupling. Consequently, the modes 11^2 and iwa are real and positive for any value of the momentum, 
and the corresponding terms in Eq. ()17p will decay exponentially in time. On the other hand, the solutions iujo and 
iui are real only if 2ga'^kFLs ~ 2TVpkpa < 1. This condition is always satisfied in the weak spin-orbit coupling 
regime, in which the diffusion equation formalism applies. We conclude that within the diffusive regime a system 
will never exhibit oscillatory dynamics. Below, we show that this conclusion holds for arbitrary spin-orbit coupling 
parameters. Nonetheless, we can formally extrapolate Eq. (|57)) outside its strict domain of validity and apply it to the 
strong spin-orbit interaction regime. In this limit, i.e. for TVFkpa > 1, the modes iujQ and iwi acquire an imaginary 
component, which generates oscillatory terms in the solution (fT7)) . We interpret this second major feature of Eq. (|57)) 
as an indication that oscillatory dynamics is a signature of strong spin orbit interactions. While the diffusion equation 
formalism cannot offer a quantitative description of this regime, it should represent a good indicator for the range of 
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FIG. 6: Time dependence of the diagonal components of the diffusion matrix in the weak spin-orbit coupling (left panel) and 
strong spin-orbit coupling (right panel) regimes. The matrix elements are evaluated for the pure Rashba case (7^ = 2, 70 = 0) 
at a momentum {k^, ky) = {0.5 /Ls, 0). In the weak spin-orbit coupling limit all the quantities decay exponentially, while in 
the opposite limit the charge and the Sy channels exhibit oscillatory behavior. In the pure Rashba (or pure Dresselhaus) case 
the out-of-plane spin component does not couple to the charge, due to symmetry. For k parallel to the x-axis, this is also the 
case for the in-plane Sx spin component (see Eq. ([Bl|). 



parameters where the oscillatory behavior is likely to appear and for the degree in which various components of the 
spin and charge sectors are affected. For example, from Eq. (|B1I) wc observe that the time evolution of the diagonal 
Sz component Dzz is controlled by the modes iu!2 and itos, which do not couple with the charge. Consequently, in the 
pure Rashba (or pure Dresselhaus) case the out-of-plane spin component does not exhibit oscillatory dynamics for any 
strength of the spin-orbit coupling. This result is connected to the isotropic nature of the spin-orbit splitting in these 
particular cases. In contrast, the dynamics of the in-plane spin components is coupled with the charge dynamics. 
Explicitly, the time dependence of the matrix element Dqq describing charge relaxation is given by 



Dooit,k) = <^ 



cosh 



sinh(i^l-4g2^2fe2 



COS 



(1^4 



g'^fifk'^ — 1 



if '^.g^ik < 1, 



if 2g^ik > 1, 



(38) 



where the wave vector k is expressed in units of 1/Ls and time in units of spin relaxation time, Tg. Notice that in the 
absence of spin-charge coupling, Dqq reduces to the standard solution of the diffusion equation, Doo{t, k) = exp(— fc^f), 
while the time dependence of the matrix elements associated with the spin degrees of freedom is determined by S2 and 
S3 only. This behavior is illustrated in Fig. [6] which shows the time dependence of all diagonal matrix elements Dait, k) 
for k = (0.5,0). The curves in the left panel were calculated in the weak spin-orbit coupling regime corresponding 
to r = 0.0005 and g ~ 1000. Notice that in the charge channel we have the standard exponential decay, exp(— fc^i), 
while in the spin channels the effect of the finite relaxation time ^ Ty ^ Tz/2 is evident. By contrast, in the 
strong-coupling limit, the charge and Sy channels exhibit oscillatory behavior. Also notice that, as a result of the 
spin-charge coupling, the charge relaxes on the same time scale as the spin degrees of freedom. We mentioned above 
that for a pure Rashba (or pure Desselhaus) spin-orbit coupling, the charge does not couple with the out-of-plane 
spin component due to symmetry. In the example shown in Fig. [6] similar symmetry reasons apply to the Sx spin 
component, due to the particular choice of momentum along the x-axis. However, in general all spin components will 
couple to the charge. For example, the general expression for the D^z matrix element describing the coupling between 
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FIG. 7: Time dependence of the diagonal components of the diffusion matrix in the weak spin-orbit coupling (left panel) and 
strong spin-orbit coupling (right panel) regimes for 7^ = 1 = —70. The matrix elements are evaluated at a momentum {kx, 
ky) = (O.5/L3, 0). In the weak spin-orbit coupling limit all the quantities decay exponentially, while in the opposite limit they 
exhibit oscillatory behavior. Notice that due to the particular choice of parameters, in the weak-coupling limit the in-plane 
spin components have identical relaxation curves, D^x = Dyy, while in the strong-coupling limit the two quantities differ due 
to the k-dependent coupling to the charge channel. 



charge and out-of-plane spin density is 

Do^{s,k) = 2^(g ~ ^^)[7i?7i3(7n/"o - 7flMi) + 2s(7r,Mi -7flMo)], (39) 

where A(s, k) is the determinant of D~^. Let us consider a set of parameters for which Dpz is non-zero. We show in 
Fig. [7] the time dependence of the diagonal matrix elements in the presence of both Rashbe and Dessselhaus terms, 
so that 7fl = 1 and 7/3 = —1. In terms of the original coupling constants, this case corresponds to a « 0.47 F, 
(3i 0.19 F, and f]^ « 0.66 F. Notice that in the strong-coupling limit (right panel in Fig. [7|) all channels exhibit 
oscillatory behavior, although the relative amplitude for the Sz component is very small. Again, we observe that Dqq 
decays on a time scale of order r^, instead of the standard exponential decay exp(— fc^t). We conclude that strong 
spin-orbit interaction can generate two effects: 1) an oscillatory structure in the decay curves, and 2) an increase in 
the decay rate of the charge channel. The first effect should be observable above a certain critical momentum of order 
kc ~ gfii/Lg. The second effect can be understood as a re- normalization of the diffusion constant. In the limit of 
small momenta equation (|38p reduces to 

I?oo(t,k) w exph(l + 2g^^lf)kh]. 

Consequently, in this limit the charge diffuses at long lengthscales according to the standard diffusion equation but 
with a re-normalized diffusion constant. This effective diffusion constant is enhanced by a factor 1 -I- 2g^iif over its 
bare value: 

«(l + 2.gV)4^. (40) 

Formally, the enhancement factor can be very large. However, we have to remember that these conclusions are based 
on an extrapolation of the diffusion equation formalism outside its domain of validity and, therefore, should be seen as 
having qualitative rather than quantitative relevance. Strictly speaking, in the strong spin-orbit coupling regime the 
diffusion equation formalism breaks down and the notion of a diffusion coefficient becomes meaningless. The coupled 
spin-orbit dynamics is no longer described by a set of differential equations of type (jlip and p3p . 
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Our calculation using the diffusion equation formalism indicates that oscillatory dynamics may appear in a system 
with strong spin-orbit interaction as a consequence of spin-charge coupling. The effective strength of this coupling 
depends not only on the overall strength of the spin-orbit interaction, i.e. on T, but also on the interplay between 
different types of spin-orbit couplings. The net result of these competing effects is contained in the expressions of the 
spin-charge coupling coefficients /io and given by Eq. (jl2[) . In the examples considered above (see Fig. [6] and [7]), 
we have (/io, fii) ~ (0, F^) for the pure Rashba case and (/io,/ii) = (1.09 F^, 1.40 F^) for the 77? = 1 = — 7d case. 
For comparison, the maximum values of the spin-charge coupling coefficients are |/io| ~ 1.10 F^ and w 1.65 F^. 
However, for other values of the spin-orbit coupling parameters these coefficients, and implicitly the spin-charge 
coupling effects, may be much weaker. In particular the special case corresponding to the symmetric points defined 
by 7i? = 7d = ±-\/2 deserves special attention (similar considerations can be made for the case 'Jr = —71? by 
interchanging the and '-' components). In this case fio ~ fii ~ and, consequently, the spin-charge coupling 
vanishes. Using the rotated coordinate system and expression (|13p for the polarizability matrix, we obtain for the 
diffuson 
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(41) 



The effects of the spin-charge coupling are totally absent at the symmetry points, regardless of the strength of the 
spin-orbit interaction. In addition, the diffusion matrix is independent of the /c_ wave-vector component and has a 
gapless mode for k+ — ±V2, as discussed in the previous section. 

To understand the effects of spin-orbit coupling in a system with arbitrary spin-orbit interaction, we calculate the 
diffusion matrix corresponding to wave- vectors along the k± directions. From Eq. psp we obtain the general expression 
for the matrix D^'^\s,k±) given in Appendix [Bl The real time behavior can be obtained by Fourier transforming 
(jB2|) with respect to the frequency. Each matrix element has a time dependence of the form given by Eq. (|17|). except 
that this time the k-vector is restricted to the k± directions. The four modes that control the coupled spin-charge 
dynamics are 
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(42) 



where /io and /ii are expressed in terms of the spin-orbit coupling parameters in Eq. \12\ . Similar expressions can 
be written for the fc_ direction, using the correspondence rules '+'^'-', 7/5 ^ID, fJ-i — /ii and iuji iLU2- As 
in the isotropic case discussed above, the charge mode iujo couples with only one spin mode, iui for the k-^- direction 
and iu;2 for the fc_ direction, while the other spin modes remain unaffected by the spin-charge coupling. Also, for 
these special k-space directions, the dynamics of the out-of-plane component of the spin is determined by the spin 
modes that do not couple to the charge, i.e. D^^ in Eq. (|B2p depends on iuj2{k+) and iw3(fc+) only. We conclude that 
oscillatory dynamics, the re-normalization of the diffusion constants, or any other spin-charge coupling effect should 
not be observable in the relaxation of an injected out-of-plane spin density with a wave-vector along the k± directions. 
However, for an arbitrary direction, the effects of the spin-charge coupling are present, as we have seen in Fig. [71 
as all the modes contribute to the spin dynamics. The set of four relaxation rate modes, has some main qualitative 
features that are independent of direction: i) there are always at least two purely real modes, one of them being iujs, 
and ii) an imaginary component may be acquired by the charge mode iujQ and one of the spin modes itoi or iw2, 
exactly which one of them depending on the direction in k-space. The appearance of the imaginary component takes 
place above a certain minimum value of the wave-vector satisfying the condition 2g\fiQ ± ^ii\kmin = 1 T 1r1d/2. We 
mention that the existence of such an imaginary component does not necessarily generates an oscillatory relaxation of 
the spin and charge densities. The general solution in this case contains two exponentially decaying terms an addition 
to the oscillatory term with exponentially decaying amplitude and the existence of an oscillatory component in the 
relaxation curves is determined by the relative intensity of these contributions. Finally, at small values of the wave- 
vector, the coupled spin-charge relaxation is always monotonic. However, the relaxation is, in general, anisotropic 
and the diffusion coefficients get re-normalized. 
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V. 



SUMMARY 



In this paper, we studied spin diffusion in a generic spin-orbit coupled system. We found that there are a number 
of different dynamic regimes and phenomena controUed by the relative values of the spin-orbit coupling parameters. 
These unusual phenomena include an enhancement of the spin polarization lifetime at finite wave- vectors, anisotropic 
spin diffusion, an oscillatory behavior of the spin relaxation in real space, an enhancement of the effective charge 
diffusion coefficient, and possible real-time oscillations in spin diffusion dynamics. The existence and manifestations 
of these phenomena depend in a non-trivial manner on the strength of the various terms contributing to the spin-orbit 
interaction. 

One of the main results of the paper is a general analytic form of the spin-charge coupled diffusion equation 
in the presence of all possible types of spin-orbit couplings. The parameters of this equation can be conveniently 
parameterized by three dimensionless coupling constants: The overall strength of the spin-orbit coupling (F) and two 
Rashba-type (7/?) and Dresselhaus-type (7d) contributions to spin precession. The last two parameters must satisfy 
the physical constraint 7^+7!) < 4. We showed that spin diffusion can be completely characterized by the pair of these 
dimensionless parameters (77^, 7d). We have identified the domains in the parameter space characterized by either 
non-zero or zero k- vector maxima of the spin relaxation spectrum, which determine the asymptotic dynamics at large 
times. One particularly important example we considered is the vicinity of the symmetry point |7/j| = |7i3| = V2, 
where the lifetime becomes arbitrarily long for certain finite k-vectors. Our analysis showed that, in order to observe 
experimentally this slow spin relaxation, special attention should be paid to minimize the cubic Dressclhaus spin-orbit 
interaction. 

We have also qualitatively considered the regime of strong spin-charge coupling and extrapolated the results of the 
diffusion equation formalism to the strong spin-orbit interaction regime. We found that the enhancement of the charge 
diffusion coefficient and the oscillatory spin relaxation observed in recent spin-grating experiments can be qualitatively 
understood as a spin-charge coupling effect. We showed that the coupling of spin and charge is in general strongly 
k-dcpendent and can be characterized by two parameters fxo and ^1, which depend on the spin-orbit interaction. 
These coupling parameters depend not only on the overall interaction strength F, but also on the interplay between 
the Rashba and Dresselhaus contributions. We argued that a necessary condition to observe an oscillatory behavior 
is a strong spin-orbit interaction, corresponding to the crossover from diffusive transport to the ballistic regime. For a 
given direction in k-space, the effect should be observable for k larger than a certain minimum wave- vector. However, 
we predict that the effect is not observable along certain directions in k-space for which spin-charge coupling is 
prevented by symmetry. For example, the out-of-plane spin polarization does not couple with the charge sector in the 
pure Rashba or pure Dresselhaus case. The coupling is also absent in general for directions in k-space corresponding 
to fc+ and k- (see text). In addition, we predict that the effect is absent if the competition between the Rashba and 
Dresselhaus terms generates small coupling constants /xo and /^i, regardless of the overall strength of the spin-orbit 
interaction. This is the case, for example, in the vicinity of the special symmetry points where /iq and /ii vanish. We 
should reiterate here that the diffusion approximation may capture spin-charge coupled dynamics only qualitatively 
if the spin-charge coupling is strong. To obtain a quantitative theory in this case, one should study the dynamics of 
the matrix distribution functions in the Boltzmann equation. 

APPENDIX A: SOLUTION OF THE SPIN DIFFUSION EQUATION IN THE LIMIT OF SMALL k- 

WAVE-NUMBERS 

The secular equation A(s,k) = has simple analytical solutions if one of the wave-vector components A:_ or 
vanishes. The asymptotic real time and real space spin dynamics is controlled by the minima of the dispersion curves 
itt'j(k) that are either at zero momentum, or at a finite value along the k^ or A;_ directions. To determine this 
asymptotic behavior, we need the solution of the secular equation in the vicinity of the minima. Assuming that fc_ is 
small, we obtain for an arbitrary value of k^ the following approximate solutions of the secular equation 



si(k) 



2 



(Al) 



2 



1 




) 



2 



+ 4kli-fn + lD)^~^2Mk+) kl+0{kt), 



where the coefficients of the quadratic contributions in fc_ are 
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\/(1-^)+4^-+(7h + 7d)^ 



Similar expressions can be obtained for wave- vectors with small and arbitrary fc_ components using the correspon- 
dence rules si 32, k+ fc_ and jd —Id- We notice that this quadratic approximation is valid away from level 
crossing points. For example, the denominator of vanishes for a set of parameters corresponding to a degeneracy 
point for the iuji and iuj2 modes. 

The quadratic expansion around the k = minimum of the iuJi mode has the form 



zc.i(k) = l-^ + A;2 +C(^)(^^„^^) kl+0{k'), (A3) 



with 



Q '{1R,1d) - 1 - -y^-y^ ■ (A4) 

2 

This zero momentum minimum exists provided Cg^^ {iri Id) > 0, i.e. for spin-orbit coupling parameters corresponding 
to the red and blue zones in Fig. [U but it docs not contribute to the dynamics along the line with equal Rashba and 
Dresselhaus couplings, = 713. Similarly, we obtain for the zero momentum minimum of the 11^2 mode that exists 
within the red and yellow regions in Fig. [T] the expression 

iu:20^)^l + ^^ + kl+cf\^R,^D) kl+0{k^), (A5) 

with 

^(2)/ N -I (7fi + 7-d)^ 

'{iR^lD) = 1 - -, ^n^o ■ (A6) 

^ 2 

For the finite momentum minimum of the iuj2 mode at A:^ that exists within the blue and green zones (see Fig. [T]) 
we obtain from equations (|Aip and (|A2p . 



i, , (Ir) ~ ^ Jl^l^ClR (1 - (7fl + 7o)^ p(2) ,0 N2 , ^(2) ,2 /A7N 

iuj2{k) ^ - + - _____ + (fc+ - + C.o (A7) 



where the coefficients of the quadratic terms are 

B^Xir.Id) = l- \ ' (A8) 
+ (7i?. + 7-d) 

^(2), ^, X -, (7fi - ID? IR + 1I + f 7fl7u - 1 
+ (7i?. + 7d)^ 7^ + 71, 4- 57/?7i5 - 1 

Similar expressions can be obtained for the finite momentum minimum of the iuji mode, corresponding to the yellow 
and green regions in Fig. [U using the correspondence rules: iuj2 iuJi, fc+ ^ fc-, and —"fD- 



APPENDIX B: EXPLICIT EXPRESSIONS OF THE DIFFUSION MATRIX IN THE PRESENCE OF 

SPIN-CHARGE COUPLING 



The diffusion matrix propagator for pure Rashba or pure Dresselhaus spin-orbit interaction can be obtained by 
inverting the polarizability matrix from Eq. (jlip . after setting = and fio = (for the Rashba case) or "fa = 
and fii = (for Dresselhaus spin-orbit interaction). We obtain 
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A(s,k) 
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A(s,k) 
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-igp,ik^ 
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k^ky{{s-l)j^ + {s+l)g^t^l) 
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(S-S2)(S-S3) 
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where 7 = ±2 for the Rashba interaction and 7 G [—2, 2] for the Dresselhaus case, the spin-charge coupHng parameter 
is given by Eq. p2p . and A(s,k) = (s — Sq){s — Si)(s — S2)(s — S3) is the determinant of the polarizabiHty matrix. 
The solutions of the secular equation A = are Sj(k) = — iwj(k) + 1 + k'^ with the dispersion relations for the modes 
itUj given by Eq. ([57)1 . 

For a general set of spin-orbit interaction parameters, one can easily invert the matrix in Eq. (|13[) and obtain the 
diffuson matrix for wave vectors oriented along the '+' and '— ' axes of the rotated coordinate system. Explicitly, for 
the fc-i- direction we have 



V 



„ TR70 
^ 2 

{s-sa}{s-si) 



{s-so){s-si) 







s+1 

(S-S2){S — S3) 



'(7H+7-D)fc + 
(s-S2)(s-S3) 



-»g(A'o+/'i)fc+ 

(s-So)(s-Si) 



(s-so)(s-si) 







{s — S2)(s-sa) 



(s-S2)(s-S3) 



(B2) 



where Sj = —i^^j + I + k'^ represent the four modes described by equation (|42|) . Notice that the charge couples only 
with the S- spin component, i.e. with the in-plane component perpendicular to the wave-vector. A similar expression 
can be obtained for the fc_ direction using the standard correspondence rules '+' <-^' — ' and jd — ^ — 7d- 
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